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Abstract. The high-frequency peaked BL Lac PG 1553+113 was observed in 2005 with the H.E.S.S. stereoscopic 
array of imaging atmospheric-Cherenkov telescopes in Namibia. Using the H.E.S.S. standard analysis, an excess 
was measured at the 4. Oct level in these observations (7.6 hours live time). Three alternative, lower-threshold 
analyses yield >5ct excesses. The observed integral flux above 200 GeV is (4.8±1.3 s t a t±1.0 S j,st) x 10 -12 cm _2 s _1 , 
and shows no evidence for variability. The measured energy spectrum is characterized by a very soft power law 
(photon index of r=4.0±0.6). Although the redshift of PG 1553+113 is unknown, there are strong indications that 
it is greater than z=0.25 and possibly larger than 2=0.78. The observed spectrum is interpreted in the context of 
VHE 7-ray absorption by the Extragalactic Background Light, and is used to place an upper limit on the redshift 
of 2<0.74. 
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1. Introduction 

PG 1553+113 was discovered in the Palomar-Green sur- 
vey of UV-excess stellar objects (Green et al. 1986) It 



is classified as a BL Lac object based on its featureless 
spectrum UMiller fc Green 19831 IFalomo fc Treves 1990|> 
and its significant (m p = 13.2-15.0) optical variability 
(Miller et al. 1988)| PG 1553+113 is well studied from 
the radio to the X-ray regime, and has been the 
subject of several multi- wavelength observation cam- 
paigns. In X-rays it has been detected by multi- 
ple observatories, with energy spectra measured by 



both BeppoSAX\(Donato et al. 2005)] and XMM Newton 
(Perlman et al. 2005) Based on its broad-band spectral 
energy distribution (SED), PG 155 3+113 is now classified 
as a high-frequency peaked BL Lac (Giommi et al. 1995)| 
similar to essentially all of the AGN detected at VHE (very 
high energy; >100 GeV) energies. 

The redshift of PG 1553+113 was initially determined 
to be z=0.36 IjMiller fc Green 1983ft . However, it was 
later claimed that this value is flawed as it was based 
on a spurious emission line, caused by a bright spot in 
the International Ultraviolet Explorer image, that was 
misidentified as Lyman-a l|Falomo fc Treves""T9 90 1 . To 
date no emission or absorption lines have been measured 
despite more than ten observation campaigns with op- 
tical instruments, including two with the 8-meter VLT 



telescopes (Wagner 2005) Estimates of the redshift can 
also be made using photometric observations. No host 
galaxy was resolved in Hubble Space Telescope (HST) im- 
ages of PG 1553+113 taken during the HST survey of 110 
BL Lac objects (Scarpa et al. 2000) Approximately 80% 
of these BL Lacs (88) have known redshifts, of which 
all 39 with z<0.25 and 2 1 of the 28 with 0.25<z <0.6 
have their hosts resolved (Sbarufatti et al. 2005) The 



HST results were recently used to set a lower limit of 



z>0.78 |(Sbarufatti et al. 2005)| for PG 1553+113. This 
redshift limit is based on the assumption of the nu- 
cleus being present in a typical host galaxy, and gen- 
eral correlation between photometrically determined red- 
shifts and those determined spectroscopically. However, 
a departure from this correlation, or the possibility of 
an atypical host galaxy is not unexpected. The possibil- 
ity of such a large redshift is of critical importance to 
VHE observations due to the absorption of VHE pho- 
tons l|Stecker et al. 19921 IHauser fc Dwek 2001|) by pair- 
production on the Extragalactic Background Light (EBL). 
This absorption, which is energy dependent and increases 
strongly with redshift, distorts the VHE energy spectra 
observed from distant objects, and may even render the 
VHE detection of very distant objects impossible. 

Although the redshift of PG 1553+113 is essentially 
unknown, it is viewed as a promising candidate for de- 



tection as a VHE emitter (Costamante & Ghisellini 2002) 
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based on its SED. However, it has not been previ- 
ously detected in the VHE regime. Only flux level up- 
per limits have been reported by the Whipple collabora- 
tion|(de la Calle Perez et al. 2003") and the Milagro group 
KWilliams 2004)| Here, evidence for >200 GeV 7-ray emis- 
sion from PG 1553+113 is presented, based on observa- 
tions made with the H.E.S.S. (High Energy Stereoscopic 
System) experiment. The data are used to set an upper 
limit on the source redshift. 

2. H.E.S.S. Detector & Observations 

The H.E.S.S. experiment is a stereoscopic system of 
four imaging atmospheric-Cherenkov telescopes located in 
Namibia. H.E.S.S. is designed to search for astrophysi- 
cal 7-ray emission above ~100 GeV, with a sensitivity 
(5a in 25 hours for a 1% Crab Nebula flux source at 
20° zenith angle) that allows for the detection of VHE 
emission from objects such as PG 1553+113 at previously 
undetect able flux levels . More details on H.E.S.S. can be 
found in |Hinton (2004) 

The H.E.S.S. observations of PG 1553+113 were made 
in May and August 2005. A total of 22 runs (~28 minutes 
each) were taken, yielding a total exposure of 9.6 hours. 
Seventeen of these runs pass selection criteria which re- 
move data for which the weather conditions were poor or 
the hardware was not functioning properly. The total ex- 
posure of the quality selected data set is 7.6 hours live 
time, and the mean zenith angle of these observations is 
40°. All these data were taken using Wobble mode, i.e. 
the source direction is positioned ±0.5° relative to the 
center of the field of view (f.o.v.) of the camera during ob- 
servations, which allows for both on-source observations 
and simultaneous estimation of the background induced 
by charged cosmic rays. 

3. Analysis Technique 

The data passing the run selection criteria are calibrated 
as detailed in Aharonian et al. (2004) and the event re- 



construction and background rejection (i.e. event selec- 
tion criteria) are performed as described in Aharonian 
et al. (2005a), with some minor improvements discussed 
in Aharonian et al. (2005b). On-source data are taken 
from a circular region of radius 9 cu t centered on the 
source. The background is estimated using events from 
a number of equal-size off-source regions offset by the 
same distance, relative to the center of the f.o.v., in 
the sky as the on-source region (for more details see 
lAharonian et al. 2005cfl . The significance of any excess 
is calculated following the method of Equation (17) in 



Li & Ma (1983) 

Table [T] shows the cuts on 
width (MRSW) and length 
(Aharonian et al. 2005a) on 2 



mean reduced scaled 
(MRSL) parameters 
(the square of the 
angular difference between the reconstructed shower 
direction and the source position), on individual image 
size (in photoclcctrons; PE), and on the distance of the 
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Table 1. The selection cuts applied to the data for the 
H.E.S.S. standard analysis and those applied for the spec- 
trum determination in this article. Cuts keeping only 
events with MRSL and MRSW greater than -2. Oct are 
used in both cases. Only images passing the distance and 
size cuts are used in the analysis, and a minimum of two 
such images are required. 



Cut 


MRSL 


MRSW 


a 2 

u cut 


Size 


Dist. 


Type 


max 


max 


max 


min 


max 




M 


M 


[deg 2 ] 


[PE] 


[mrad] 


Standard 


2.0 


0.9 


0.0125 


80 


35 


Spectrum 


1.3 


0.9 


0.02 


40 


35 



Table 2. Results of the H.E.S.S. observations of 
PG 1553+113 using both the standard analysis cuts and 
the lower-threshold spectral analysis cuts. Shown are the 
post-cuts energy threshold at the mean zenith angle of ob- 
servations (40°), the number of on-source and off-source 
events passing the cuts, the normalization for the off- 
source events, the observed excess (A) from PG 1553+113 
and the significance (S) of the excess. The significance re- 
sulting from the spectral analysis cuts is subject to a small 
trials factor. Therefore only the significance resulting from 
the H.E.S.S. standard analysis should be considered. 



Cut 


E t h 










S 


Type 


[GeV] 


On 


Off 


Norm 


A 


M 


Standard 


310 


457 


4126 


0.0907 


83 


4.0 


Spectrum 


230 


1807 


12753 


0.1236 


230 


5.3 



image center of gravity from the center of the f.o.v., used 
in this analysis. The standard H.E.S.S. analysis cuts 
shown in Table ^ are optimized a priori using Monte 
Carlo 7-ray simulations and unrelated off-source data 
to yield the maximum expected significance per hour of 
observation for a weak (10% Crab flux source with photon 
index T=2.6) source. Generally, the significance expected 
(and observed) is not strongly dependent on the exact 
values of the cuts. However, the size cut increases the 
energy threshold of the analysis, therefore lower values 
are more appropriate for a steep spectrum source. 



4. Results 

Table|21shows the results from the analysis of the H.E.S.S. 
observations of PG 1553+113. Using the standard cuts, an 
excess of 83 events (4. Oct) is observed from the direction 
of PG 1553+113. Figure [I] shows the on-source and nor- 
malized off-source distributions of 9 2 for the H.E.S.S. ob- 
servations. The background is approximately flat in 9 2 as 
expected, and there is an excess at small values of 9 2 cor- 
responding to the observed signal. The shape of this excess 
is consistent with that expected from a point source given 
the H.E.S.S. point-spread function. 




0.15 

e?[deg 2 ] 

Fig. 1. The distribution of 9 2 for on-source events (points) 
and normalized off-source events (shaded) from observa- 
tions of PG 1553+113. The dashed line represents the cut 
on 9 2 applied to the data. 
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Fig. 2. The energy spectrum of PG 1553+113. The hori- 
zontal error bars indicate the energy bin size. The dashed 
line represents the best \ 2 fit of a power law. 



Assuming a power-law spectrum (dN/dE ~ E~ T ) 
with a photon index of T=4.0 (justified later), the 
observed integral flux above 200 GeV is I(>200 
GeV)=(4.8±1.3 stat ±1.0 S2/s t)xlO- 12 cnr 2 s _1 . This corre- 
sponds to -2% of I(>200 GeV) determined by H.E.S.S. 
for the Crab Nebula. The flux is well below all previously 
published upper limits for PG 1553+113. No evidence for 
variability is found as a fit to the nightly integral flux 
versus time is consistent with being constant (x 2 =3.7, 7 
degrees of freedom). Additionally, quick- look results pro- 
vided by the RXTE/ASM team show no indication of flar- 
ing behavior at X-ray energies (2 to 10 keV), suggesting 
that PG 1553+113 was not in a high-emission state. 

Given the low number of excess events found using 
the H.E.S.S. standard analysis, the determination of an 
energy spectrum is not possible. Therefore, a set of lower- 
threshold cuts were generated to increase the selection effi- 
ciency for low-energy 7-rays, allowing for a spectral deter- 
mination. The cuts were generated using the exact same 
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5. Discussion 
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Fig. 3. The intrinsic photon index of PG 1553+113, with 
lcr statistical error contours, versus redshift for the cases of 
minimal and maximal EBL absorption. An intrinsic pho- 
ton index below the horizontal line (Pi n t=1.5) is not con- 
sidered realistic. The uppermost (dashed-dotted) curve is 
the sum of the intrinsic photon index and 2a statistical 
error for the case of minimal EBL absorption. 



procedure as for the standard analysis cuts, differing only 
in that they were optimized for a 1% Crab flux (>100 
GeV) source with photon index T=5.0. Using these cuts 
(also shown in TableQJ an excess of 230 events (5.3cx) from 
the direction of PG 1553+113 is found. The significance of 
this excess is clearly subject to a small trials factor and 
is therefore reduced. To be conservative, only the signifi- 
cance of the standard analysis should be considered. The 
energy spectrum determined from the data set, using the 
low-threshold analysis, is shown in Figure The best \ 2 
fit of a power law to these data yields a photon index 
r=4.0±0.6 stat , and a x 2 01 1-8 for 2 degrees of freedom. 
The systematic error on the photon index is small com- 
pared to the statistical error. 

It should be noted that the integral flux (again assum- 
ing r=4.0) determined with the "spectrum" cuts, I(>200 
GeV) = (5.8±l.l stat ±1.2 syst )xl(r 12 cm -2 s -1 , is consis- 
tent with the flux obtained using the standard cuts. This 
indicates that the increased excess found with these cuts 
is indeed expected and not the result of over-optimization 
or statistical fluctuations. In addition, a cross check was 
performed on the same data as published by H.E.S.S. on 



PKS 2005-489 |(Aharonian et al. 2005b)| Here an analy- 
sis using the spectral-determination cuts results in the 
same spectrum (r=4.0) as previously published, as well 
as the increase of the detected excess by similar fac- 
tor as for PG 1553+113. Finally, an alternative anal- 
ysis using an independent calibration and two differ- 
ent event reconstruction methods Ijde Naurois et al. "2 003 
|Lemoine-Goumard fc D cgrang e~2004| ) verifies the ob- 
served signal (>5i7 excess) and yields compatible spectra. 



PG 1553+113 possibly represents the most distant BL 
Lac object detected at VHE energies. Therefore, the ob- 
served VHE spectrum is expected to be affected by EBL 
absorption. Unfortunately the unknown redshift of this 
object does not allow its VHE spectrum to be used in 
constraining the poorly-measured EBL flux density (see 
e.g. lAharonian et al.l 2005d). However, assuming a red- 
shift allows the determination of the spectrum intrinsic 
to the BL Lac (i.e. with the effect of EBL absorption re- 
moved). The intrinsic flux is related to the observed flux by 
Fint = e T ( z ' E ' F bs, where r(z, E) is the optical depth due 



to pair production (7 VHE 7 EBL — > e + e~). Assuming that 
the main uncertainty arises from the overall normalization 
of the EBL density, this optical depth can be expressed 
as t(z, E) — f t (z, E), where r (z, E) is calculated using 
a reference EBL-density model, and / is a scaling fac- 
tor for this density. At energies relevant to the H.E.S.S. 
observations, t(z, E) increases with both redshift and en- 
ergy. Therefore, a larger redshift, or higher EBL flux den- 
sity, causes a stronger absorption of higher-energy 7-rays, 
resulting in a harder intrinsic spectrum (F; nt oc E~ Tini ). 
The EBL density model of Primack et al. (2001) is cho- 
sen here since it includes the effects of galaxy evolu- 
tion which become important at redshifts greater than 
~0.2. This model's fluxes are above recently published 
upper limits IjAharonian et all 2005d). However, scaling 
its values by /=0.85 approximately traces these limits 1 
and is therefore considered a maximal EBL flux density. 
Given that the EBL flux density is poorly-known, a min- 
imal EBL scenario is also considered. Here, the origi- 
nal model is scaled by /=0.6, approximately tracing the 
lower limits on the EBL flux density provided by galaxy 
counts (Madau & Pozzetti 2000) Using these parameter- 



izations, the intrinsic spectrum of PG 1553+113 is calcu- 
lated for different redshift values. Figure 01 shows the in- 
trinsic photon index (r; nt ) versus redshift. Assuming the 
intrinsic photon index of the BL Lac is not harder than 
Pi nt =1.5 fsee lAharonian et al.l 2005d for arguments sup- 
porting this assumption) , an upper limit on the redshift of 
PG 1553+113 of z<0.74 is determined using the minimal 
EBL. Above this redshift, where ri nt =0.3±0.6, the sum 
of Tint and twice its statistical error is harder than 1.5. 
Clearly, the redshift upper limit is smaller if a higher EBL 
flux density is assumed. This is the first time that a VHE 
observation is used to usefully constrain the redshift of an 
astrophysical object. Future observations of PG 1553+113 
should allow for a more precise determination of the pho- 
ton index, possibly resulting in a stronger redshift con- 
straint. 

Acknowledgements. The support of the Namibian authorities 
and of the University of Namibia in facilitating the construc- 
tion and operation of H.E.S.S. is gratefully acknowledged, 
as is the support by the German Ministry for Education 



1 A different model and hence scaling factor is used in 
lAharonian et al.l (2005d) . 



F. Aharonian et al.: Evidence for VHE 7-ray Emission from the Distant BL Lac PG 1553+113 



5 



and Research (BMBF), the Max-Planck-Society, the French 
Ministry for Research, the CNRS-IN2P3 and the Astroparticle 
Interdisciplinary Programme of the CNRS, the U.K. Particle 
Physics and Astronomy Research Council (PPARC), the IPNP 
of the Charles University, the South African Department of 
Science and Technology and National Research Foundation, 
and by the University of Namibia. We appreciate the excel- 
lent work of the technical support staff in Berlin, Durham, 
Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia 
in the construction and operation of the equipment. 

References 

Aharonian, F., et al. (H.E.S.S. Collaboration), 2004, Astropart 
Phys, 22, 109 

Aharonian, F., et al. (H.E.S.S. Collaboration), 2005a, A&A, 
430, 865 

Aharonian, F., et al. (H.E.S.S. Collaboration), 2005b, A&A, 
436, L17 

Aharonian, F., et al. (H.E.S.S. Collaboration), 2005c, A&A, 
441, 465 

Aharonian, F., et al. (H.E.S.S. Collaboration), 2005d, submit- 
ted for publication arXiv:astro-ph/0508073 
Costamante, L. & Ghisellini C, 2002, A&A, 384, 56 
de la Calle Perez, I., Bond, I.H., Boyle, P.J., et al., 2003, ApJ, 
599, 909 

de Naurois, M., Guy, J., Djannati-Atai, A., et al., 2003, Proc. 

of the 28 th ICRC (Tsukuba), 2907 
Donato, D., Sambruna, R.M., & Gliozzi, M., 2005, A&A, 433, 

1163 

Falomo, R. & Treves, A., 1990, PASP, 102, 1120 

Giommi, P., Ansari, S.G., & Micol, A., 1995, A&AS, 109, 267 

Green, R.F., Schmidt, M., & Liebert, J., 1986, ApJS, 61, 305 

Hauser, M.G. & Dwek, E., 2001, ARA&A, 39, 249 

Hinton, J., 2004, New Astron Rev, 48, 331 

Lemoine-Goumard, M. & Degrange, B., 2004, AIP Conference 

Proceedings, 745, 697 
Li, T. & Ma, Y., 1983, ApJ, 272, 317 

Miller, H.R., Carini, M.T., Gaston, B.J., et al., 1988, Proc. of 

the IUE Symposium (Greenbelt), 303 
Miller, H.R. & Green, R.F., 1983, BAAS, 15, 957 
Perlman, E.S., Madejski, G., Georganopoulos, M., et al., 2005, 

ApJ, 625, 727 
Madau, P. & Pozzetti, L., 2000, MNRAS, 312, L9 
Primack, J.R., Somerville, R.S., Bullock, J.S., et al., 2001, AIP 

Conference Proceedings, 558, 463 
Sbarufatti, B., Treves, A., & Falomo, R., 2005, ApJ, 635, 173 
Scarpa, R., Urry, CM., Falomo, R., et al., 2000, ApJ, 532, 740 
Stecker, F.W., de Jager, O.C., & Salamon, M.H., 1992, ApJ, 

390, 49 

Wagner, S., 2005, Private Communication 

Williams, D., 2004, AIP Conference Proceedings, 745, 499 



